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We investigate the photoproduction of p mesons in ultraperipheral heavy ion collisions at RHIC 
and LHC energies in the dipole approach and within two phenomenological models based on the the 
Color Glass Condensate (CGC) formalism. We estimate the integrated cross section and rapidity 
distribution for meson production and compare our predictions with the data from the STAR collab- 
oration. In particular, we demonstrate that the total cross section at RHIC is strongly dependent on 
the energy behavior of the dipole-target cross section at low energies, which is not well determined 
in the dipole approach. In contrast, the predictions at midrapidities at RHIC and in the full rapidity 
at LHC are under theoretical control and can be used to test the QCD dynamics at high energies. 
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I. INTRODUCTION 

The Large Hadron Collider (LHC) at CERN will start 
up in few months. Currently, there is a great expectation 
that LHC will discover the Higgs boson and whatever 
new physics beyond the Standard Model that may ac- 
company it, such as supersymmetry or extra dimensions 
However, several questions remain open in the Stan- 
dard Model, which will be probed in a new kinematical 
regime at the LHC and will determine the background for 
new physics. In particular, the description of the high- 
energy regime of the Quantum Chromodynamics (QCD) 
still is a subject of intense debate (For recent reviews see 
Ref. 1). 

In recent years we have proposed the analysis of co- 
herent collisions in hadronic interactions as an alterna- 
tive way to study the QCD dynamics at high energies 
@, H i, i, 0, S, H . The basic idea in coherent hadronic 
collisions is that the total cross section for a given pro- 
cess can be factorized in terms of the equivalent flux 
of photons into the hadron projectile and the photon- 
photon or photon-target production cross section. The 
main advantage of using colliding hadrons and nuclear 
beams for studying photon induced interactions is the 
high equivalent photon energies and luminosities that can 
be achieved at existing and future accelerators (For a re- 
view see Ref. [lCj). Consequently, studies of -fp interac- 
tions at the LHC could provide valuable information on 
the QCD dynamics at high energies. 

The photon-hadron interactions in hadron-hadron col- 
lisions can be divided into exclusive and inclusive re- 
actions. In the first case, a certain particle is pro- 
duced while the target remains in the ground state (or 
is only internally excited). On the other hand, in inclu- 
sive interactions the particle produced is accompanied 
by one or more particles from the breakup of the tar- 
get. The typical examples of these processes are the ex- 



clusive vector meson production and the inclusive heavy 
quark production, described by the processes 7/1 — > Vh 
{V = p, J/*,T) and jh -> XY (X = cc.bb), respec- 
tively. Recently, we have discussed both processes con- 
sidering pp H, H, Q , pA Q and AA d, Q collisions as an 
alternative to constrain the QCD d yna mics at high ener- 
gies (For recent reviews see Refs. dHHH). Our results 
demonstrate that their detection is feasible at the LHC. 
However, some of these results were obtained consider- 
ing the saturation model proposed by Golec-Biernat and 
Wusthoff (GBW) several years ago [l3j]. That model has 
been improved recently considering the current state-of- 
the-art of saturation physics: the Color Glass Condensate 
(CGC) formalism [14|, [l5|, [l6| . This fact motivates a revi- 
sion of some our previous estimates. In particular, in this 
paper we revise our predictions for p-meson photoproduc- 
tion in coherent AA collisions. We show that the rapidity 
distribution at RHIC is strongly dependent on the behav- 
ior of dipole-nucleus cross section at the threshold energy 
(and, therefore, model dependent), whereas the theoreti- 
cal predictions based on CGC physics at LHC energy are 
under control. Moreover, we demonstrate that the STAR 
Collaboration data for the rapidity distribution can be 
described by a phenomenological model which captures 
the basic properties of the quantum evolution in the CGC 
formalism. 



This paper is organized as follows. In the next section 
we present a brief review of coherent hadron-hadron in- 
teractions, introducing the main formulas and discussing 
the QCD dynamics at high energies. In particular, we 
present the saturation models considered in the calcula- 
tions. In Section IIIII we compare our results with the ex- 
perimental data from STAR Collaboration at RHIC and 
present our predictions for the p production at LHC. Fi- 
nally, in Section IIVI we summarize our main results and 
conclusions. 



II. COHERENT AA INTERACTIONS AND QCD 
DYNAMICS AT HIGH ENERGIES 



Lets consider the hadron-hadron interaction at large 
impact parameter (b > 2Ra) at ultra relativistic energies. 
In this regime we expect the electromagnetic interaction 
to be dominant. In heavy ion colliders, the heavy nuclei 
give rise to strong electromagnetic fields due to the coher- 
ent action of all protons in the nucleus, which can interact 
with each other. The photon emmited from the electro- 
magnetic field of one of the two colliding hadrons can in- 
teract with one photon of the other hadron (two-photon 
process) or can interact directly with the other hadron 
(photon-hadron process). In what follows our main focus 
will be in photon - hadron processes. Considering the re- 
quirement that photoproduction is not accompanied by 
hadronic interaction (ultra-peripheral collision) an ana- 
lytic approximation for the equivalent photon flux of a 
nuclei can be calculated, which is given by [To| 



dN 1 (u) _ 2 Z 2 a e 



dw 



(1) 



where ui is the photon energy, 7^ is the Lorentz boost of a 
single beam and Ko(fj) and K\(ff) are the modified Bessel 
functions. Moreover, fj = 2uu Ra/jl a,ndU(fj) = K\ (fj) — 
Kq (fj) . The Eq. fT} will be used in our calculations of p 
meson production in A A collisions at RHIC and LHC. 

The cross section for the p meson photoproduction in 
a coherent hadron-hadron collision is given by, 



a(AA -> AAp°) = 2 dw dt 



dN^(w) da 
dw dt 



(W 1P ,t) 12) 



where 4nr is the differential cross section for the process 



dt 

(jA -> p°A), w mi „ = M 2 /4 7L m p , W 2 p = 2wy/S^m 
\J Snn is the c.m.s energy of the hadron-hadron system. 
The coherence condition limits the photon virtuality to 
very low values, which implies that for most purposes, 
they can be considered as real. In particular, for PbPb 
collisions at LHC the Lorentz factor is jl — 2930, giving 
the maximum c.m.s. jN energy W lp ~ 950 GeV. The 
factor two in Eq. {2]| takes into account the fact that 
nuclei can act as both target and photon emitter. The 
experimental separation for such events is relatively easy, 
as photon emission is coherent over the entire nucleus 
and the photon is colorless we expect the events to be 
characterized by intact recoiled nuclei (tagged nuclei) and 
a two rapidity gaps pattern (For a detailed discussion see 

We describe the vector meson production in the color 
dipole frame, in which most of the energy is carried by the 
hadron, while the photon has just enough energy to disso- 
ciate into a quark-antiquark pair before the scattering. In 
this representation the probing projectile fluctuates into 
a quark-antiquark pair (a dipole) with transverse separa- 
tion r long after the interaction, which then scatters off 
the hadron [l7| . In the dipole picture the imaginary part 
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FIG. 1: (Color online) Predictions for the rapidity distribu- 
tion of p° photoproduction at RHIC energy considering dis- 
tinct theoretical approaches. Data from STAR Collaboration 



of am plit ude for vector meson production reads as (see 
Refs. MEGl) 

A(x,A)=f dzd^be-^ 1 -^ ^^,)^^) 

where ^ ^(z, r) and ty^-Az, r) are the light-cone wave 
function of the photon and of the vector meson, respec- 
tively, with the quark and antiquark helicities labeled by 
h and h. The variable r defines the relative transverse 
separation of the pair (dipole), z (1 — z) is the longitu- 
dinal momentum fractions of the quark (antiquark), A 
denotes the transverse momentum lost by the outgoing 
proton (t = —A 2 ), and x is the Bjorken variable. The 
transverse distance from the center of the target to one of 
the qq pair of the dipole is denoted by b. The differential 
dipole-target cross section is referred as -^g- The pho- 
ton wave functions appearing in Eq. §S§ are well known 
in literature [l9| . For the meson wave function, we con- 
sider the Gauss-LC model [l9| (the results are insensitive 
to a different choice). Moreover, we assume an effective 
light quark mass (m u — 0.14 GeV). Finally, the differ- 
ential cross section for vector meson photoproduction is 
given by 



da 
~dt 



1 

16^ 



|^(x,A)| 2 (l +/ 3 2 ), (4) 



where f3 is the ratio of real to imaginary parts of the 
scattering amplitude. For the case of rho meson photo- 
production, skewness corrections are not important and 
they are not taken into account. A recent discussion on 
the vector meson production in eA colliders and a de- 
tailed list of references about this subject can be found 
in Ref. 

The differential dipole-target cross section contains all 
information about the target and the strong interaction 
physics. In the Color Glass Condensate (CGC) formal- 
ism [l4|, [H, EH j h encodes all the information about the 
non-linear and quantum effects in the hadron wave func- 



3 



tion. It can be obtained by solving an appropriate evo- 
lution equation in the rapidity y = ln(l/a;) and its main 
properties are: (a) for the interaction of a small dipole 
(r -C 1/Qs) the dipole-target satisfies the color trans- 
parency property (cr q q cx r 2 ), which characterizes that 
the system is weakly interacting; (b) for a large dipole 
(r 3> 1/Qs), the system is strongly absorbed. This prop- 
erty is associate to the large density of saturated glu- 
ons in the hadron wave function. In our analysis we 
will consider two saturation models generalized for the 
photon-nucleus interaction. The first one, generalizes the 
Iancu-Itakura-Munier (IIM) model [21| . introducing the 
impact parameter dependence. In this model the differ- 
ential dipole-nucleus cross section is given by [HJ 



IIM 



d 2 b 



= 2\ I - exp 



-~AT A (b)a™(x,r) 



where T^(b) is the nuclear profile function (obtained from 
a 3-parameter Fermi distribution for the nuclear density); 
the dipole-nucleon cross section is parametrized as [21| 



2 n 2 \ 7eff (x,r) 



*ZM=*a\^^r'''' i for rQ s < 2 , 
1 — exp [— c In (drQ s )\ , for rQ s > 2 , 



where the saturation scale is given by Q s {x) = {xo/x) x ^ 2 , 
with Ao = 0.7. The constants c and d are determined 
from the condition that <7di P and its derivative with re- 
spect to rQ s are continuous at rQ s — 2. The expression 
for rQ s > 2 (saturation region) has the correct func- 
tional form, as obtained from the theory of the Color 
Glass Condensate (CGC) [14[. Moreover, -f c ff (x, r) = 
7sat + la ^/xy 3 ^ i s the effective anomalous dimension, 
which determines the behavior of the dipole cross section 
in the color transparency regime (rQ s < 2) and intro- 
duces evolution effects not included in the GBW model 
[13]. The parameters of model are determined from a fit 
to i*2 data. The IIM model can be considered a phe- 
nomenological model for the quantum limit of the CGC 
formalism, since it captures the basic properties of the 
quantum evolution, describing both the bremsstrahhmg 
limit of linear small- a; evolution (BFKL equation) as well 
as nonlinear renormalization group at high parton densi- 
ties (very small- a;). For comparison we also consider the 
IP-SAT model generalized for nuclei target proposed in 
Ref. [23j, which is a a phenomenological model for the 
classical limit of the CGC. In this model, the differential 
dipole-nucleus cross section takes the form 



d?b 



= 2 



I — exp 



2N, 



■a.(ji 2 )xg{x,ii 3 )AT(b) 



where the scale p 2 is related to the dipole size r by 
p 2 = 4/r 2 + /j,q and the gluon density is evolved from 
a scale /Zq up to p 2 using LO DGLAP evolution without 
quarks, assuming that the initial gluon density is given 
by xg(x, fify = A g x~ x q (1 — a;) 5 - 6 . The values of the pa- 
rameters /ig, A g and X g are determined from a fit to F% 
data. The parameter set used in our calculations is the 



one presented in the first line of Table III of [l9(. It is 
important to emphasize that both models describe quite 
well the small- a; HERA data (x < 0.01). Therefore, the 
study of other obscrvablcs which are strongly dependent 
on is very important to constrain the underlying 
QCD dynamics at high energies. 

An important comment is in order here. The descrip- 
tion of the vector meson production using the color dipole 
picture is justified at small values of x (x < 0.01) and 
large x contributions are not included. As we will show 
in next section, in the kinematical region of large rapidi- 
ties at RHIC, the dipole-cross section is probed in the 
large- a; region. In order to take this fact into account, 
we will multiply the differential dipole-nucleus cross sec- 
tions presented above by a threshold factor, (1 — x) n , as 
performed in previous studies [1, [l|[ . Notice that in the 
IP-SAT model, a large- a; correction is already included in 
the proton gluon distribution function. We discuss the 
consequence of this phenomenological procedure in what 
follows. 



III. RESULTS 

Lets calculate the rapidity distribution and total cross 
section for the p° meson photoproduction in coherent 
A A collisions. The distribution on rapidity y of the 
produced final state can be directly computed from Eq. 
(|2|), by using its relation with the photon energy to, i.e. 
y oc In (2u}/Mp). Explicitly, the rapidity distribution is 
written down as, 



dN^juj) 
: ^ — dij — a i A ~*P° A v w ) W 



da [A + A -> A (g) p° <g> A] 
dy 

where ® represents the presence of a rapidity gap. Con- 
sequently, given the photon flux, the rapidity distribution 
is thus a direct measure of the photoproduction cross sec- 
tion for a given energy. 

Initially, lets investigate the p° photoproduction at 
RHIC energies. The present calculation update our pre- 
vious calculation 0] and allow us to compare our pre- 
dictions with those from Refs. 1251 . |26( and the experi- 
mental data presented in Ref. [2J] . In Fig. Q] we com- 
pare our predictions with those from Klein and Nystrand 
(KN) [1| and Frankfurt, Strikman and Zhalov (FSZ) 
[26j |. which are based on distinct theoretical approaches 
(For a detailed discussion of these models see [5]). Wc 
verify that the IP-SAT model overestimates the STAR 
data at mid rapidity, whereas the IIM model seems to 
be consistent with the overall normalization of the ex- 
perimental measurements. In comparison with the KN 
and FSZ predictions, the predictions of the saturation 
models are not strongly suppressed at mid rapidities and 
■s/s = 200 GeV. Such suppression occurs only at smaller 
energies. 

The results presented in Fig. [1] were obtained using 
n = 4 in the threshold factor. In order to verify the 
dependence of our predictions in the model used for the 



4 



(a) IPSAT 




>- 

"D 
73 



(b) MM 



/// 



STAR 
n = 3.0 
n = 4.0 
n = 5.0 
n = 8.0 



\\\ 



FIG. 2: (Color online) The rapidity distribution of p° meson photoproduction at RHIC energy (y^i = 200 GeV) for IP-SAT 
(left panel) and IIM (right panel) models. Data from STAR Collaboration [24j]- 
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FIG. 3: (Color online) Dependence on the threshold factor of 
the rapidity distribution of p° photoproduction at the LHC 
(PbPb collisions and y/s = 5.5 TeV) considering the IIM 
model. 



low energy regime, which is not under theoretical control 
in the dipole approach, we show in detail in Fig. [2] the 
dependence on the threshold factor, (1 — x) n , taking dif- 
ferent values for the index n (n = 3,4,5,8). The form 
of the threshold factor was inspired by the quark count- 
ing rules. The effect is completely negligible at central 
rapidity, whereas it is important at forward rapidities. 
The dependence on n starts to be very important for ra- 
pidities \y\ > 2. It implies that our predictions in the full 
rapidity range are strongly dependent on the model used 
to describe the low energy regime. This fact has impor- 
tant consequences when computing the total exclusive 
cross section at RHIC energy. In contrast, we have that 
at LHC energies the threshold factor is negligible for most 
of rapidity range, being of marginal importance for very 
forward rapidities, as can be observed in Fig. [3j where 
we present the IIM predictions for the rapidity distribu- 
tion for p° photoproduction in PbPb collisions at LHC 
(y^s = 5500 GeV). A similar behavior is observed using 



the IP-SAT model. Therefore, the color dipole model is 
a suitable approach to describe p meson production at 
LHC. 

In Fig. 0] we present our predictions for the rapid- 
ity distribution for the LHC energy considering PbPb 
and CaCa collisions. The results for both models are 
shown for the full rapidity range. The rapidity distri- 
bution at LHC probes a large interval of photon-nucleus 
center of mass energy since W^ A ~ M p ^fse v , which 
corresponds to very small x ~ M p e~ y / ^fs. Therefore, 
its experimental analyzes can be useful to determine the 
QCD dynamics. As already observed at RHIC energies, 
the IIM predictions are smaller than the IP-SAT one for 
collisions with both nuclei. However, the difference be- 
tween the predictions of the models increases with the 
energy, as observed in the left panel of Fig. [5] where we 
present the energy dependence of the total cross section 
calculated assuming the cut |y| < 1. The reason to im- 
pose this kinematical cut is based on the sensitivity to 
the large- x correction to the dipole cross section. The 
plot is motivated by recent studies presenting the en- 
ergy dependence of the p° photoproduction cross section 
for the RHIC energies. The large difference between the 
predictions is directly associated to the distinct energy 
dependence predicted by the IIM and IP-SAT model for 
the photo-nucleus cross section. In Fig. [5] (right panel) 
we present our predictions for the energy dependence of 
da/dy at central rapidity (y = 0), which is not depen- 
dent on the model used to describe the large y region. A 
similar behavior is also observed for this quantity, which 
implies that its experimental study can be useful to dis- 
criminate between the saturation models. Our predic- 
tions for the rapidity distribution at y — and LHC 
energies (See F ig. jU) can be compared with those pre- 
sented in Refs. [25|, |2g|. In comparison with Ref. (25j. our 
prediction for p° photoproduction in CaCa collisions is a 
factor > 3 larger, which is directly associated to the dis- 
tinct energy dependence among the models. In Ref. [25[ 



5 



J3 

> 

o 



. IIM(Ca) 

LrlU IPSAT(Ca) 

IIM(Pb) 

IPSAT(Pb) 





-8 -7 -6 -5 -4 -3 



-I I 

Y 



3 4 5 6 7 g 



FIG. 4: (Color online) Predictions for the rapidity distribu- 
tion of p° photoproduction in PbPb and CaCa collisions at the 
LHC. The results for IP-SAT and IIM models are presented. 



A 


IIM 


IP-SAT 


Ca 


395 mb (16985.0) 


572 mb (24600.0) 


Pb 


5977 mb (2510.0) 


8596 mb (3610.0) 



TABLE I: The integrated cross section (events rate/second) 
for the p° photoproduction in AA collisions at LHC energies 
and different nuclei (A — Pb and Ca). 



it is assumed that the p° photoproduction is described by 
soft physics, which implies a smooth energy dependence. 
In contrast, in the models based on CGC physics, it is 
expected that at high energies the saturation scale de- 
termines the QCD dynamics implying a stronger energy 
dependence. In comparison with Ref. (2d| . our predic- 
tions for da/dy(y = 0) for p° photoproduction in PbPb 
collisions are similar. 

We are now ready to compute the total exclusive cross 
section, Eq. @, and production rates for LHC energies 
considering the distinct phenomenological models. The 
results are presented in Table [H At LHC we assume the 
design luminosity £ L hc = 0.42 (43.0) mb^s" 1 for PbPb 
(CaCa) collisions. The corresponding cross sections are 
very large, which implies that the experimental analyzes 
of this process at LHC should be feasible. In particular, 
as the p° meson photoproduction is an exclusive reaction, 
the separation of the signal from hadronic background 
would be relatively clear. Namely, the characteristic fea- 
tures are low px p° spectra and a double rapidity gap pat- 
tern. Moreover, the detection (Roman pots) of the scat- 
tered nuclei can be an additional useful feature. More- 
over, in comparison with the hadroproduction of vector 
mesons, the event multiplicity for photoproduction inter- 
actions is lower, which implies that it may be used as a 
separation factor between these processes (For a detailed 
discussion see Refs. HHi]). In 

comparison with previous 



results 0, l25l \2dl1 our predictions are smaller than those 
presented in [26| . obtained using the Generalized Vector 
Dominance Model (GVDM), with the IIM one being a 
factor « 1.6 smaller. In contrast, they are a factor > 3 
larger than those presented in [25j . where a VDM plus 
a classical mechanical calculation for nuclear scattering 
was used in order to calculate the rapidity distribution 
and cross sections (For a detailed comparison of these 
models with the dipole approach see Section IV in Ref. 
5]). In comparison with our previous results [j| the IIM 
and IP-SAT predictions for PbPb collisions are smaller 
than those obtained using the GBW model, which is di- 
rectly associated to the distinct energy dependence of the 
dipole-target cross sections. 



IV. CONCLUSIONS 

The QCD dynamics at high energies is of utmost im- 
portance for building a realistic description of AA col- 
lisions at LHC. In this limit QCD evolution leads to 
a system with high gluon density, denoted Color Glass 
Condensate (CGC). If such system there exists at high 
energies it can be proven in coherent AA collisions at the 
LHC. In this paper we have analyzed the p° meson pho- 
toproduction considering the dipole approach and most 
recent phenomenological models based in CGC physics. 
We shown that the predictions of the dipole approach 
for rapidity distribution at RHIC are strongly dependent 
on the behavior of dipole-nucleus cross section at thresh- 
old energy (and, therefore, model dependent). It implies 
that this approach is not able to give reliable predictions 
for the total cross section at RHIC. However, at central 
rapidities (y < |1|), where are the experimental data of 
the STAR Collaboration, our prediction are not strongly 
dependent on the model used at low energies. At this 
kinematical region we demonstrate that the STAR data 
can be described by the IIM model, which captures the 
basic properties of the quantum evolution in the CGC 
formalism. Furthermore, we demonstrate that the pre- 
dictions for the rapidity distribution and total cross sec- 
tion based on CGC physics are under theoretical control 
at the LHC energy. Moreover, we predict very large cross 
sections and event rates, which implies the experimental 
study of this process is feasible at LHC. Consequently, 
we believe that it can be used to study the QCD at high 
energies. 
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FIG. 5: (Color online) The total exclusive cross section as a function of energy in the range 50 < y/s < 500 GeV (left panel) 
for the IP-SAT and IIM models (using cut \y\ < 1). The differential cross section da/dy at central rapidity y — as a function 
of energy (right panel) . 
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